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ABSTRACT 

We have surveyed six distinct 'empty fields' using the Austraha Telescope Compact Ar- 
ray in an ultra-compact configuration with the aim of imaging, with a high brightness 
sensitivity, any arcmin-scale brightness-temperature anisotropies in the background 
radio sky. The six well-separated regions were observed at a fi-equency of 8.7 GHz 
and the survey regions were limited by the ATCA primary beams which have a full 
width at half maximum of 6 arcmin at this frequency; all fields were observed with a 
resolution of 2 arcmin and an rms thermal noise of 24 /iJy beam~^. After subtracting 
foreground confusion detected in higher resolution images of the fields, residual fluctu- 
ations in Stokes I images are consistent with the expectations from thermal noise and 
weaker (unidentified) foreground sources; the Stokes Q and U images are consistent 
with expectations from thermal noise. 

Within the sensitivity of our observations, we have no reason to believe that there 
are any Sunyaev-Zeldovich holes in the microwave sky surveyed. Assuming Gaussian- 
form CMB anisotropy with a 'flat' spectrum, we derive 95 per cent confidence upper 
limits of Q flat < 10-11 /iK in polarized intensity and Q fiat < 25 /xK in total intensity. 
The ATCA filter function peaks at I — 4700 and has half maximum values at I = 3350 
and 6050. 
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1 INTRODUCTION 

Today, we are fairly confident that the large scale struc- 
ture grew out of primordial seed perturbations by gravita- 
tional instabilites in a dark matter dominated universe. The 
physics that dominates the coupling between matter inho- 
mogeneities and the cosmic microwave background (CMB) 
radiation anisotropies is different on different angular scales 
and the multifarious coupling mechanisms are believed to 
result in a complicated spectrum of radiation anisotropies 
(see, for example, Hu, Sugiyama & Silk 1996 for a recent 
review) . 

The primary anisotropies in the CMB that result from 
sub-horizon couplings via gravitational, Doppler and acous- 
tic effects are expected to create features in the anisotropy 
spectrum on angular scales exceeding several arcmin; the 
cutoff at small angular scales depends on photon diffusion, 
recombination physics, any reionization at late epochs and 
the cosmological parameters that determine the angular 
size-redshift relation. Improvements in the precision of mea- 
surements of the band powers in anisotropies in the CMB, 
combined with large-scale-structure observations, lead to re- 
finements in the parameters of the background cosmology, 
models for structure generation and matter content of the 
universe (see Bond & Jaffe 1997, Lawrence, Scott & White 
1998 and Gawiser & Silk 1998 for examples of the values 



of CMB anisotropy measurements in constraining structure 
formation models). 

The power in the CMB anisotropy spectrum towards 
small angular scales is expected to be strongly damped: the 
spectral cutoff may constrain the ionization history of the 
primordial gas and mechanisms for reionization (Dodelson & 
Jubas 1995). Reionization may erase primary anisotropies on 
small angular scale; however, via the 'Vishniac effect', reion- 
ization may generate secondary fluctuations of signiflcant 
magnitude (Hu, Scott & Silk 1994). Gravitational lensing 
by large-scale structure may also transfer anisotropy power 
from larger to smaller scales enhancing the amplitude of pri- 
mary anisotropies at small scales (Metcalf & Silk 1997) . The 
time dependence of the gravitational potential during late 
non-linear structure formation may imprint anisotopies on 
the CMB via the 'Rees-Sciama effect' and this may be the 
dominant anisotropy on sub-arcminute scales (Seljak 1996). 
Another plausible source for small-angular-scale anisotropy 
is the Sunyaev-Zeldovich effect in the relativistic plasma co- 
coons of radio galaxies (Yamada, Sugiyama & Silk 1999). 
Anisotropies in the CMB on arcmin scales could also pro- 
vide a powerful constraint on any cosmological magnetic 
field (Subramanian & Barrow 1998). 

If structure forms by the growth of initial small per- 
turbations, the CMB anisotropies are expected to be par- 
tially polarized. This is because the matter inhomogeneities 
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cause the radiation to be anisotropic, and Thomson scatter- 
ing of the anisotropic radiation at the last scattering sur- 
face generates hnear polarization (Bond & Efstathiou 1984; 
Hu &: White 1997). The fractional polarization in the pri- 
mary anisotropies is expected to be greater at smaller an- 
gular scales and the detection of polarization in the CMB 
anisotropies may be a fundamental test of, and provide an 
additional observational constraint on theories for structure 
formation. 

The secondary anisotropies result from interactions in 
the post-rccombination era and during non-linear struc- 
ture formation. An important effect here is the Sunyaev- 
Zeldovich (S-Z) scattering of the CMB as it encounters hot 
gas in the deep potential wells of clusters of galaxies (see 
Birkinshaw 1998 for a recent review). The expectations for 
the CMB anisotropy arising from this effect are difficult to 
compute, and depend on the abundances of dark matter 
halos with cluster masses that are formed at different cos- 
mic epochs, and the processes by which the gas accumu- 
lates in these potential wells and is heated. Nevertheless, 
using models for the evolution of the intracluster medium 
and estimates of the abundances of cluster potentials from 
the Press-Schechter formalism or self-similar evolutionary 
models, distribution functions have been derived for the S- 
Z sources and these have been used to estimate the CMB 
anisotropies in different cosmological models and scenarios 
for the evolution of the gas (Bartlett & Silk 1994; Marke- 
vitch et al. 1994). Blind surveys for S-Z clusters are expected 
to provide valuable inputs to our understanding of the late 
stages of structure formation. 

A deep radio survey with the VLA in an HST medium 
deep survey field does appear to have detected a radio decre- 
ment (Richards et al. 1997); independently, a Ryle-telescope 
observation towards a pair of quasars has also shown a 'hole' 
in the radio sky (Jones et al. 1997). In both cases, follow-up 
in the optical and X-ray bands have not revealed any cluster 
in these directions. Because the Compton y parameter which 
determines the magnitude of any S-Z decrement towards a 
hot-gas cluster is distance independent, the S-Z flux den- 
sity (the integrated decrement) may be observable to higher 
redshifts as compared to the optical/X-ray luminosities. The 
absence of luminous counterparts in other bands indicates 
that any clusters that may be responsible for the decrements 
via the Sunyaev-Zeldovich effect should be at moderately 
high redshifts: the existence of collapsed massive structures 
at early epochs defies conventional wisdom (Bartlett, Blan- 
chard & Barbosa 1998) and has led to alternate hypotheses 
for microwave decrements (Natarajan & Sigurdsson 1999). 

In 1990, we commenced a programme to search for 
arcmin-scale anisotropies with the Australia Telescope Com- 
pact Array (referred to as the ATCA experiment). The first 
results based on a deep survey of a single field were reported 
in Subrahmanyan et al. (1993); the sensitivity was improved 
with further observations using improved receivers and the 
derived upper limits were reported in Subrahmanyan et al. 
(1998). We have now used the Australia Telescope to im- 
prove the sky area covered and imaged six independent sky 
regions with high brightness temperature sensitivity. Our 
goal was to make a blind search for any 'holes' in the mi- 
crowave sky apart from deriving estimates of the arcmin 
scale CMB anisotropy in total intensity and polarization. In 
this paper, we present these observations. 



2 ATCA OBSERVATIONS 

The Australia Telescope Compact Array (ATCA; see The 
Australia Telescope 1992) is an E-W array with five 22-m 
diameter antennae movable on a railtrack. The best surface 
brightness sensitivity is achieved when the array is in a close- 
packed configuration and we have configured the ATCA in 
a 122-m array, with the five antennae equispaced along the 
E-W track on stations 30.6 m apart, for these observations. 
In this configuration, observations at declinations lower than 
—50° would be shadowed in some hour angles and shadowing 
results in spurious correlations in the baseline pairs involv- 
ing antennae whose projected areas overlap. The brightness 
sensitivity in the observations would progressively decrease 
for fields at higher declinations because of the increase in 
projected baselines at large hour angles. The 22-m antennae 
cannot be placed on stations closer than 30.6 m; therefore, 
we have made all observations in fields located about —50° 
in declination so that the projected spacings between adja- 
cent antennae would be close to the shadowing limit of 22 m 
over a large range of hour angles. The observations would 
have the best brightness sensitivity without shadowing re- 
lated problems. 

Most physical mechanisms for primary anisotropics in 
the CMB result in wavelength-independent fluctuations in 
the CMB thermodynamic temperature. At centimetre wave- 
lengths, where the sky brightness of the 2.73 K CMB radiar- 
tion may be approximated by the Rayleigh- Jean's law, these 
CMB brightness temperature fluctuations may be consid- 
ered to be wavelength independent. On small angular scales 
the dominant cause of CMB temperature anisotropy may be 
the S-Z effect: this has a wavelength dependence. However, 
for this effect, a fairly constant temperature decrement is 
expected at centimetre wavelengths. 

The anisotropy power in the CMB sky is usually de- 
scribed in terms of C;'s in multipole Z-space, where Ci = {\ 
Cim P) and Cim's are the coefficients in the spherical har- 
monic expansion of the fractional temperature anisotropy 
on the celestial sphere. The anisotropy variance, which is ex- 
pected in any differential measurement of the sky brightness 
temperature, and is made with an i-space window function 
of octave width, will be about l{l + l)Ci; therefore, the CMB 
anisotropy spectrum is usually a plot of l{l + l)Ci versus 
As discussed earlier, the anisotropy power due to primary 
anisotropies is expected to cut off at small angles and con- 
sequently on arcmin scales, the expected variance in most 
experimental efforts may be expected to be lower if the effec- 
tive beam is smaller. There have been estimates of the CMB 
anisotropy spectrum arising from the decrements due to the 
thermal S-Z effect in a cosmological distribution of clusters 
containing hot gas; the resulting l{l + l)Ci is fairly flat in 
the /-space regime corresponding to arcmin scales (Molnar 
1998). 

In any image made with the ATCA in its ultra com- 
pact array configuration with the five 22-m antennae located 
30.6 m apart along the E-W track, provided the observations 
are made with same times and using receivers with the same 
system temperatures, the brightness sensitivity (i.e., the rms 
noise in the image in units of brightness temperature) would 
be independent of observing frequency. Observations made 
at higher frequencies would have smaller beams and the win- 
dow function corresponding to the measurements would be 
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located higher up in Z-space. Consequently, ATCA observa- 
tions made at higher frequencies would be attempting to 
make measurements of the CMB spectrum at larger but 
would be less likely to detect any primary anisotropics. On 
the other hand, ATCA observations at all frequencies would 
be equally sensitive to 'flat band' anisotropics. 

The ATCA currently operates in the 20, 13, 6 and 3 cm 
wavelength bands; in its ultra-compact configuration, the 
images at all those wavelengths would be confusion lim- 
ited at levels above those required to detect most CMB 
anisotropics. The confusion at these wavelengths arises from 
unresolved extragalactic synchrotron sources which typically 
have a steep flux-density spectrum with a mean spectral in- 
dex a ~ —0.6 define the spectral index by the relation 
5*1^ ~ v"). These sources have an integral number count that 
approximately follows the relation iV(> S) ~ v~°'^ S~^ '^Q., 
where S is the flux density, v is the observing frequency 
and O is the beam solid angle. We may assume that the 
discrete sources are Poisson random distributed on the sky 
and that the rms confusion noise is roughly the flux den- 
sity at which we may expect one source in every synthesized 
beam. Because the synthesized beam scales as f2 ~ , we 
may expect the rms confusion noise to scale as AS ~ . 
The brightness temperature scales as T ~ S/{Q,v^) and 
it follows that the rms brightness temperature 'noise' on 
the image would scale as AT ~ . Every octave in- 
crease in frequency reduces the confusion noise by a factor 
about 4.6; therefore, we have opted to make the ATCA CMB 
anisotropy observations in the highest available 3 cm band. 

When observing with the ATCA in the ultra compact 
122-m array configuration, with five antennae equi-spaced 
30.6-m apart along the E-W railtrack, images made using 
just the four 30.6-m baselines — instead of using all ten 
baselines — give the highest surface brightness sensitivity. 
Images made using just the four 30.6-m baselines between 
adjacent a ntenn ae would have an rms thermal noise that is 
a factor 10/4 larger, and a synthetic beam that is at least 
a factor four larger in solid angle. Because the brightness 
temperature sensitivity depends on the rms thermal noise 
(in units of flux density per beam) and inversely on the beam 
solid angle, the images made omitting baselines exceeding 
30.6 m would have a brightness sensitivity that is at least 
a factor 4/-^10/4 better. For this reason, the images that 
represent the CMB anisotropy search have been constructed 
using the four 30.6-m baselines alone; images made using the 
remaining six baselines exceeding 30.6 m have been used to 
derive the foreground discrete source confusion in the fields. 

As discussed in section 6.1 below, the CMB anisotropy 
search using the four 30.6-m baselines is sensitive to 
anisotropy power in an /-space window defined by the tele- 
scope filter function. The observations are sensitive to CMB 
anisotropy with spatial frequencies that are covered by the 
30.6-m baseline and, therefore, the filter function would have 
a smaller fractional bandwidth (and appear to be narrower) 
in the case of observations closer to the celestial pole where 
the projected baseline is, on the average, larger. As stated 
above, the observations are restricted to declinations south 
of —50'' to avoid shadowing. Consequently, the brightness 
sensitivity of the 30.6-m interferometer, as also the frac- 
tional bandwidth covered by the interferometer in /-space, 
is a maximum for observations at declination —50°. For ob- 
servations at declination —50°, the 30.6-m baseline would 



have a projected baseline close to the dish diameter of 22 m 
for a significant fraction of the observing time; consequently, 
the filter function would not be narrow. 

The half power width of the synthesized beam for these 
low-resolution images is about 2 arcmin for 3-cm observa- 
tions at —50° declination. This beam size is more or less 
right for the detection of S-Z anisotropies due to clusters at 
moderate redshifts. 

2.1 Field selection criteria 

Four regions were initially selected at about —50° decli- 
nation, where highest-surface-brightness observations were 
possible in the 122-m array without shadowing. These were 
located well off the Galactic plane. The Molonglo 408 MHz 
survey (Large et al. 1981) lists sources with flux densities 
exceeding 0.7 Jy at this declination and is believed to be 
complete to 1 Jy, the PMN survey at 4.85 GHz (Wright et 
al. 1994) lists sources exceeding 40 mjy with 90 per cent 
completeness: the selected regions had no listed sources in 
either of these catalogues within 30 arcmin of their centres. 
There are also no strong sources listed in either of these 
catalogues that are in the vicinity of the selected regions. 
We have examined the DSS data of the UK Schmidt IllaJ 
plates of these regions to verify that there are no obvious 
peculiarities in the optical — no obvious nebulosities, clus- 
ters, or higher than usual star counts. The coordinates of 
the centres of the selected regions are in Table 1. Two pairs 
of regions were selected that were about 12 hr apart in right 
ascension, so that the deep surveys for CMB anisotropies 
could be made in 12-hr LST ranges that did not overlap, 
and the observations could use the entire time the ATCA 
was specially configured in the 122-m array for the CMB 
anisotropy experiment. 

These four regions were first observed at 20 cm with 
the ATCA in its 750 m array configuration. The aim was 
to detect weaker confusing sources in the regions for the 
purpose of using the images to select a set of small fields 
for deep 3 cm imaging for CMB anisotropy. The observa- 
tions were made during November 14-17, 1996, in the 0.75A 
configuration and using two 64-MHz wide bands centred at 
1344 and 1434 MHz. Visibility measurements over spatial 
wavelengths 360-3400 A may be expected to be sensitive to 
confusing sources with angular scales up to a few arcmin and 
these images would be useful in selecting the fields for CMB 
anisotropy searches on arcmin scales. The absolute flux den- 
sity scale was set by short observations of PKS 1934—638 
whose flux density was adopted to be 15.01 and 14.83 Jy re- 
spectively at 1344 and 1434 MHz. Each region was observed 
for about 11 hr using dual linear polarizations; the arcmin 
resolution images had rms noise of about 0.17 mJy beam"'^ 
before any correction was made for the attenuation due to 
the primary beam. 

Sources may be reliably detected in the images with 
flux density exceeding about 0.5 mJy beam~^. The ATCA 
primary beam has half power width of about 33 arcmin at 
21 cm. Adopting an integral source count 

iV(> S) = 1966S'7y'' ™ sr"^ (1) 

at 1.4 GHz, based on the differential source counts in the 
flux density range 1-50 mJy compiled in Windhorst et al. 
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Table 1. Regions surveyed for confusion- free fields. 

Region RA(J2000) DEC(J2000) Peak Flux density Total Flux density 

at 20 cm at 20 cm 

A 03 28 26.40 -.51 00 00.00 16.6 mjy bcam^i 145 mjy 

B 01 43 26.40 -50 48 21.60 19.5 mjy beam"! 205 mjy 

C 14 39 33.84 -50 10 51.42 16.6 mjy beam"! 199 mjy 

D 11 45 23.04 -50 57 08.64 21.1 mjy beam"! 207 mjy 

The peak and total flux densities at 20 cm are quoted for images not corrected 
for the attenuation due to the primary beam. 



(1985), regions with 33 arcmin diameter may be expected to 
have about 2 sources exceeding 25 rn,Jy. The peak and total 
flux densities in the 20 cm images of the individual regions 
are also listed in Table 1. No surveyed region has sources 
exceeding 25 mJy implying that the four regions surveyed 
have relatively moderate numbers of the stronger sources. 

The imaging for CMB anisotropies was done in the 3 cm 
wavelength band where the ATCA primary beam drops to 
10 per cent of its peak value at a radial distance of 5 arcmin. 
In any region with radius 5 arcmin, about 3 sources may be 
expected on the average at 21 cm wavelength with flux den- 
sity exceeding 0.5 mJy. We identified field centres within the 
four regions surveyed that do not have sources with 21 cm 
flux density exceeding about 0.5 mJy within 5 arcmin radii. 
Assuming a mean spectral index a = —0.5, the fields may 
not be expected to have sources exceeding about 0.2 mJy 
at 8.7 GHz. A short ATCA observation at 8.7 GHz revealed 
a > 0.5 mJy source in a potential sub-field in region A — 
an inverted spectrum source — and this field centre was 
consequently rejected. 

A total of six field centres were selected for the deep 
survey at 8.7 GHz, their coordinate centres are listed in Tar- 
ble 2. 

2.2 Observing strategy 

The strategy adopted was the same as that in previous 
ATCA searches for CMB anisotropy (Subrahmamyaii et al. 
1993; Subrahmanyan et al. 1998). The observations were 
made during July 10-29, 1997, in the 122B configuration 
of the ATCA. The five ATCA antennae provided four 31- 
m baselines, three 61-m baselines, two 92-m baselines and 
one 122-m baseline. All fields were observed in full Earth- 
rotation Fourier-synthesis mode using the E-W array in 12- 
hr sessions so as to obtain visibilities along complete ellip- 
tical («, t;)-tracks in the spatial frequency domain. Observa- 
tions were made in a pair of 128-MHz wide bands centred 
at 8640 and 8768 MHz. All four polarization products — 
XX, YY, XY and YX — were measured for every antenna 
pair between the X and Y orthogonal linearly polarized sig- 
nals from each antenna. The visibility data were recorded in 
spectral line mode with 32 channels covering the 128 MHz 
bands. 

The flux density scale was set by short observations of 
PKS 1934—638 made daily; the flux density of this primary 
calibrator was adopted to be 2.84 & 2.79 Jy respectively 
at 8640 and 8768 MHz. Secondary calibrator sources were 
observed in short 3-min scans every 20 min during the ob- 
servations. The secondary calibrators and the flelds were 
both observed with the antennae pointing at the calibra- 



Table 2. Fields observed for CMB anisotropies. 



Field 


RA(J2000) 


DEC(J2000) 


Observing Time 


A2 


03 27 07.68 


-50 56 22.0 


44 hr 


Bl 


01 43 50.84 


-50 57 55.0 


36 hr 


B2 


01 42 04.21 


-51 00 47.0 


47 hr 


CI 


14 39 36.72 


-50 15 15.0 


45 hr 


Dl 


11 46 07.71 


-51 00 55.0 


46 hr 


D2 


11 44 38.95 


-50 55 37.0 


46 hr 



tor sources/field centres, but with the array phased towards 
a sky position offset 2"^ in right ascension and 30 arcmin 
in declination. This offset between the pointing and phas- 
ing centres is not expected to cause any significant visibility 
amplitude de-correlation within the spectral channels; on 
the other hand, systematic errors of various kinds that of- 
ten appear at the phase centre in Fourier-synthesis images 
would, in our case, be about 36 arcmin from the field cen- 
tre (pointing centre) and several primary beamwidths away. 
The visibility data were time averaged 30 s online. 

2.3 Imaging 

The visibility data were imaged using the miriad data re- 
duction software. 

Because the ATCA correlator has an inherent triangular 
weighting for the lag-spectrum in the 128-MHz bandwidth 
mode, alternate channels axe redundant in the 33-channel 
frequency spectrum. Any offsets in the signal sampler levels, 
or coupling of the sampler clock birdies into the RF signal, 
are expected to lead to errors in the end channels of the 
spectrum, and these may spill-over to the even numbered 
channels. For these reasons, we rejected all even numbered 
channels and the end channels. We also rejected channels 
numbered 3 and 31 because the signal-to-noise ratio was 
poor in these channels close to the band edges. A total of 
13 useful spectral channels were retained giving an effective 
bandwidth of 104 MHz. 

The phase centres were changed off-line to coincide with 
the antenna pointing centres (field centres) by appropriate 
corrections to the channel visibility phase for the secondary 
calibrator and field data. The complex antenna gains, the 
complex spectral response (bandpass) and the polarization 
leakages were derived from the secondary calibrator obser- 
vations and these were used to calibrate the field visibil- 
ities. The field visibilities were then averaged over the 13 
frequency channels and transformed into the four Stokes par- 
rameters. 

Images were made from concatenated data bases sepa- 
rately for each of the six fields. No tapering was applied to 
the visibility data and the 30-s visibilities were weighted 
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Figure 1. 
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of the six fields in Stokes I made using visibilities in the 0-1 kA range. Contours at 24.2 /ijy bcam^-'^ X (—10, —8, 
, 10, 12, 14, 16). In this figure, as also in all the following figures, the locus of the first null in the primary beam 
dotted circle with a 7 arcmin radius. 
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based on the recorded system temperatures to minimize 
the image noise. Hereinafter, we refer to images made us- 
ing the 30.6-m basehne visibilities (with spatial wavelengths 
less than 1 kA) and rejecting all spacings other than those 
between adjacently located antennae as LR images. Those 
images made using exclusively the longer spacings exceeding 
1 kA are referred to as HR images. All images were synthe- 
sized with 10 arcsec pixels and 512 pixels a side, the choice 
ensured that the synthesized beams wore well ovorsamplcd 
and the total area imaged well exceeded the primary beam 
size. While synthesizing each LR and HR image, the cor- 
responding synthetic beam pattern images were also made 
and these point-spread-function images are referred to as 
LR-PSF and HR-PSF images. The images were not decon- 
volved. 

In Table 2 are listed the total observing time on each of 
the six fields, the harmonic mean time per field was 43.6 hr. 
The average system temperature of the 3 cm ATCA system 
at the observing frequency was about 42 K (Gough 1994). 
Given that the observations were made in dual polarization 
and in a pair of bands each with an effective bandwidth of 
104 MHz, and adopting values of 0.88 for the correlator effi- 
ciency and 0.64 for the aperture efficiency of the ATCA an- 
tennae, it follows that the 'naturally weighted' images made 
using all the 10 baselines would be expected to have an rms 
noise of about 15.1 /itJy beam~^. The LR imag es wo uld be 
expected to have an rms noise that is a factor 10/4 larger, 
i.e., 23.9 /iJy bcam^^, and the HR images would be ex- 
pected to have an rms noise that is a factor ^10/6 larger, 
i.e., 19.5 jiiy beam~^. 

LR images of the six fields, made in Stokes I, are dis- 
played in Fig. 1 as contour images. HR images of these fields 
made in Stokes I are displayed in Fig 2. The images have 
not been corrected for the attenuation due to the primary 
beam and have not been deconvolved. The primary beam 
regions are indicated in the images as dotted circles repre- 
senting the location of the first null of the primary beams, 
these have radii of 7 arcmin. We have derived the weighted 
rms intensity in these images by weighting the image pixel 
intensities by the value of the primary beam at the pixel 
locations. These image rms computed for the individual LR 
images and HR images are listed in Table. 3. 

Most radio sources in the sky have extremely small, if 
any, circular polarization. Therefore, Stokes V images are 
expected to be almost entirely thermal noise images. This 
noise is uniform over synthesized images and does not follow 
any taper due to the primary beam. We have constructed 
Stokes-V LR and HR images of the six fields and computed 
the image rms noise over large areas of the individual im- 
ages. The rms noise values are listed in Table 3. These 
represent our best estimates of the thermal noise compo- 
nent of the jHxcl variances in the images in Fig. 1 and 2. 
It may be noted that the mean variance in these Stokes 
V images is (24.2 /itJy beam~^)^ for the LR images and 
(19.5 ^Jy beam"^)'^ for the HR images. Those arc within 3 
per cent of our estimates made above of the thermal noise 
based on our knowledge of the telescope parameters. 

3 FOREGROUND CONFUSION 

Clearly, the image rms in Stokes I exceeds that expected 
from thermal noise (the Stokes V image rms values); at least 



in the case of the fields A2, Bl, B2 and Dl some part of the 
excess variance will be due to the discrete sources seen in 
the corresponding HR images. 

3.1 Subtracting confusion models 

We have developed an algorithm for modelling the discrete 
sources in the fields based on the HR images and used these 
to subtract a significant part of the confusion. The inputs 
to the algorithm are the LR and HR images of any field as 
well as the corresponding synthetic beams: the LR-PSF and 
HR-PSF images. 

The algorithm iteratively identifies point-source com- 
ponents in residual HR images. Whenever point sources are 
to be subtracted from any image, the corresponding syn- 
thesized beam (PSF) image, scaled in amplitude to have a 
peak equal to the point-source flux density and shifted in sky 
position to the location of the component, is subtracted. 

In each iteration, first, a residual HR image is made 
by subtracting the component list identified in the previous 
iteration. Second, estimates are derived for the rms intensity 
in source-free regions using pixels more than 14 arcmin from 
the centre. No further components are added to the model 
and the iterations are stopped if 

(i) the number of components in the model is six, or 

(ii) the peak in the residual HR image ^ four times the 
rms in source free areas, or 

(iii) the image rms within 14 arcmin of the centre is within 
10 per cent of the rms in source-free regions and no peak 
exceeds 5 times this rms. 

New components arc added to the component list in each 
iteration by locating the peak pixel in the residual HR im- 
age within 20 arcmin of the centre. An initial estimate of the 
fiux density and position of this new component is derived by 
fitting a 7 X 7-pixel HR-PSF beam patch to a 7 x 7-pixel im- 
age patch around the position of the peak. Having updated 
the component list with this estimate for the parameters of 
the new component, the iteration ends by optimizing all the 
component parameters jointly by minimizing the variance 
in a residual HR image obtained by subtracting the compo- 
nents. 

In Table 4, we have listed the components identified 
in the six fields. It may be noted here that all components 
were constrained to be positive, and the estimation of com- 
ponent parameters did not use the LR images in any way. 
The residual LR and HR images obtained by subtracting 
these components are shown in Figs 3 & 4. 

The rms of the pixel intensities in the residual LR and 
HR images, computed after weighting the pixels by the pri- 
mary beam, are also listed in Table 3. It may be noted that 
the variance in these residual images are, in every case, in 
excess of the variance in the corresponding Stokes V images. 
The field Bl appears to have a significant negative feature 
close to the image centre. Visual comparisons between the 
residual LR and HR residual images do not show any corre- 
spondence in the positions of the peaks in the images. 

The variance in the Stokes V images are a measure 
of the thermal noise in the images of the various fields. 
We have computed the weighted mean of the variances in 
the Stokes I images of the six fields, separately for the 
LR and HR images, using the Stokes V image variances 
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Figure 2. HR images of tlie six fields in Stokes I made using visibilities exceeding 1 kA. Contours at 19.5 /iJy beam ^ x (—6, —4, —2, 
2, 4, 6, 8, 10, 12, 14, 16, 18, 20). 
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Table 3. Rms pixel intensities in the observed fields. 









Residual 






Residual 


Field 


LR image 


LR image 


LR image 


HR image 


HR image 


HR image 




Stokes I 


Stokes V 


Stokes I 


Stokes I 


Stokes V 


Stokes I 


A2 


97.4 


23.9 


52.3 


47.9 


19.1 


27.8 


Bl 


121.8 


26.9 


87.9 


49.6 


21.3 


28.5 


B2 


152.9 


23.0 


41.1 


56.5 


18.3 


20.9 


CI 


48.3 


25.5 


48.3 


25.2 


20.4 


25.2 


Dl 


165.0 


23.0 


46.0 


80.3 


19.1 


27.3 


D2 


54.7 


22.8 


30.3 


24.3 


18.5 


20.1 



All rms noise values are in /ijy beam ^ . The rms in Stokes I images were computed 
from primary-beam weighted intensities. 



Table 4. Confusion models for the observed fields. 



Field 


Component flux density 


RA offset 


DEC offset 




(MJy) 


(arcmin) 


(arcmin) 


A2 


183 


-0.4 


-0.3 




169 


1.4 


-1.2 




107 


7.8 


1.8 


Bl 


253 


1.1 


-0.4 


B2 


321 


0.4 


0.9 




163 


-0.6 


-0.3 




99 


-8.3 


-3.9 




90 


-9.1 


-16.9 


01 








Dl 


455 


-0.8 


0.2 


D2 


94 


-0.4 


1.8 



No sources were subtracted from field CI. 

as weights. Before subtracting the confusion models, the 
mean variance was (51.7 /jJy beam"'^)^ in the HR and 
(117.1 yuJy beam-i)2 for the LR images; after confusion 
subtraction, the mean variances were (25.0 /iJy boam"^)'^ 
for the HR and (52.0 /iJy bcam^^)'^ for the LR images. The 
weighted mean thermal noise (Stokes V image variance) is 
(19.4 /iJy beam-i)^ for the HR and (24.1 /xJy beam-^)^ for 
the LR images. The 'excess' variance is (15.8 /ttJy beam~^)^ 
for the HR and (46.1 /uJy beam~^)^ for the LR images. 

3.2 Residual confusion 

Windhorst et al. (1993) have derived the radio source counts 
at /ijy flux densities using deep images made at 8.44 GHz 
with the VLA. Adopting a mean spectral index a = —0.5 at 
these frequencies, the differential source counts are 

^ = (-5.05 ± 0.80) X 10-^^4-o%^jf -' /"Jy"' axcmin-^ (2) 

when scaled to 8.7 GHz. 540 ;jjy denotes the flux density in 
units of 40 pJy and the errors in the normalization and in 
the index are for a power-law fit that has the flux density S 
normalized to 40 pJy. 

We have simulated skies with a Poisson random dis- 
tribution of faint sources adopting the source counts given 
in equation 2 with the nominal values quoted for the nor- 
malization and index. Areas of the sky of about 2 arcmin 
in diameter, roughly corresponding to the half-power areas 
of the synthesized beams of the LR images, would be ex- 
pected to have one source on the average exceeding 70 /i.]y. 
In the simulations, sources in the flux density range 0.1- 
500 /uJy were included: because the differential source counts 



scale as (dN/dS) - 5'-^■^ the variance contribution from 
sources in logarithmic intervals in flux density would scale 
as cr^ ~ 3° '^ , implying that weak sources below this range 
— which may be expected to be present in large numbers 
in every synthesized beam — would contribute little to the 
image variance. Our fields are selected to be devoid of strong 
sources exceeding 0.5 mjy. 

For the simulations, weighted-mean synthesized beam 
patterns were computed by combining all the visibilities 
(of all the six fields) and computing LR and HR synthe- 
sized beams for the combined data base. In each simulation, 
the sources were convolved by these beam patterns to pro- 
duce separate LR and HR images. Thermal noise, whose 
autocorrelation was consistent with the square of the cor- 
responding synthesized beam patterns, and with variance 
(19.4 /jjy beam"^)^ for the HR and (24.1 /iJy beam"^)^ for 
the LR images, were added to the simulated images. 

The simulated images were processed by the same algo- 
rithm, described in section 3.1, to subtract discrete source 
confusion. The residual LR images had a mean variance of 
(58.8 fj,3j beam"^)^ whereas the HR images had a mean 
variance of (24.5 pJy beam"^)'^. We have grouped the sim- 
ulations ill lots of six, to simulate our observations of six 
independent fields, and computed the mean residual vari- 
ance in the groups of six fields. Their distribution implies 
that the observed rms pixel intensity, when computed as 
an average over six residual LR fields, may be expected to 
have a ±l-a spread of 50.8-65.8 ^Jy beam^^. The corre- 
sponding ±l-cr spread in the HR image rms values is 23.2- 
25.6 /uJy beam~^. The observed values of the rms pixel in- 
tensity as computed from the residual LR and HR images 
are 52.0 /xJy beam^^ and 25.0 pJy beam^^ and lie within 
the expected il-cr spreads. 

The conclusion we make from the simulations is that 
the 'excess' variance we observe in the LR as well as HR 
images, after subtracting the derived confusion models, is 
consistent with the 'excess' expected from undetected weak 
sources. All of the 'excess' variance may be attributed to 
unsubtracted confusion. 

4 FIELD Bl 

The residual image of field Bl, in Fig. 3, shows a sig- 
nificant negative close to the field centre. A fit to the LR 
visibility data yields a value of —202 /iJy and a position 
offset —17" in RA and —28" in Dec from the field centre. 
The peak of the negative feature is a factor 7.5 times the 
rms thermal noise. The corresponding residual HR image of 
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Figure 3. LR images of the six fields in Stokes I with the sources listed in Table 4 subtracted. Contours at 24.2 /xjy beam ^ x (—8, 
—7, —6, —5, —4, —3, —2, 2, 3, 4, 5). Using equation 10, we may regard these images as representing contours of brightness temperature 
at 28.7 AtK x (-8, -7, -6, -5, -4, -3, -2, 2, 3, 4, 5). 
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Figure 4. HR images of the six fields in Stokes I after subtracting the sources listed in Table 4. Contours at 19.5 /iJy beam~^ x (—4, 
—3, —2, 2, 3, 4). From the telescope filter function corresponding to the HR beam, we derive the relation AS = Q.lbQ fi^t- Using this 
factor, the contours may be regarded as being at 131 /^K X (—4, —3, —2, 2, 3, 4). 
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Figure 5. LR and HR images of the Bl_OFFSET field. The field centre is shifted 2.2 arcmin from the centre of the Bl field. No sources 
have been subtracted. Contours at 40 /ijy bcam^^ X (—6, —5, —4, —3, —2, 2, 3, 4, 5, 6). The position of the possible negative decrement 
is indicated by the star symbol in the LR image. We may regard these images as representing contours of brightness temperature at (—6, 
-5, -4, -3, -2, 2, 3, 4, 5, 6) x 47 fiK for the LR image and 270 nK for the HR image. 



Field Bl, shown in Fig. 4, shows no evidence for any negar- 

tive at this sky position: the HR image value at the position 
of the negative (in the LR image) is 70 ^iJy beam^^. 

It may be noted here that the negative peak in the LR 
image of the residual Bl field (Fig. 3) does not appear to 
have the sidelobe pattern expected for an unresolved nega- 
tive source: this means that either the negative source has 
structure on this scale or that the negative arises from a 
chance superposition of negative sidelobes of weak confus- 
ing sources in the field. 

With the aim of examining the reality of this fea- 
ture, we made a separate observation of the field with the 
ATCA antennae pointed towards position RA(J2000) = 
Q^h 43m g^^gg^ Dec.(J2000) = -50° 58' 27'.'33. This field 
centre was offset 2!2 (-127" in RA and -32" in Dec) from 
the previous pointing centre. We refer to this field centre 
as Bl.OFFSET. The 253 ^Jy source detected in the HR 
image of the Bl field is now 3.2 arcmin from the field cen- 
tre and may be expected to appear in the HR image of the 
Bl.OFFSET field with a flux density of 120 ^Jy. Because 
the shift corresponds to about 40 per cent of the primary 
beam FWHM, we may expect a significant change in the im- 
age fluctuations that arise from the numerous weak sources 
in the sky. The position of the negative peak, which was seen 
in the residual Bl image, should now be 1.8 arcmin from the 
new field centre. Correcting for the attenuation due to the 
primary beam, we may expect the negative to appear with 
a value —159 /uJy beam~^. 

We display HR and LR images of the Bl .OFFSET field 
in Fig. 5. No sources have been subtracted from these im- 
ages. The rms thermal noise is 53 /iJy beam^^ in the LR 
image and 42 /iJy beam^^ in the HR image as estimated 
from the corresponding Stokes V images. At the position in 
the HR image, where we expect a 120 /iJy source, the pixel 
intensity is 73 /iJy beam~^. At the position of the possible 
decrement (indicated by the star symbol in Fig. 5), the pixel 
intensity is —45 //Jy beam~^. The negative is not seen at the 



2(7 level. 

We have examined the residual LR images created by 
our simulations to infer the decrements that may be ex- 
pected from the negative sidelobes of residual unsubtracted 
positive confusing sources in the fields. We find that one in 
six fields may be expected to have a decrement exceeding 
—250 /iJy beam~^. A decrement of the magnitude observed 
in Field Bl, which is —206 /iJy beam^^, may be expected in 
35 per cent of the observed fields. The simulations indicate 
that the observed decrement in Bl may well be owing to 
residual confusion. 

Because the Bl_OFFSET image has not confirmed the 
decrement, and because the decrement in the Bl field does 
not have the sidelobe structure expected, we are led to ques- 
tion the reality of the decrement. We hypothesize that the 
decrement probably arises as a result of (i) weak confusion 
in the LR Bl image, and (ii) as a result of errors in the pa- 
rameter estimates of the strong confusing source seen in the 
HR image of field Bl: these errors may be due to confusion 
in the HR image. 



5 HOLES IN THE MICROWAVE SKY? 

In any sky image made at foHz GHz with a beam which has 
a half-maximum size of 6 arcmin x <p arcmin, the measured 

flux density decrement (ASjy) at any image pixel may be 
converted to brightness temperature units using the relation: 



AT = 4.6 



8.6 ' 


2 


1 


fcHz 




ex (j) 



In Table 5, we have listed the peak decrement detected in 
each of the six fields observed. The values are listed for the 
residual LR images and are computed over the half-power 
areas of the primary beams (regions of 3 arcmin radius) and 
using images that have been corrected for the attenuation 
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Table 5. Peak decrements ob- 



served 


in the residual 


LR im- 


ages. 






Field 


(/iJy beam~l) 


{^lK) 


A2 


-178 


-254 


Bl 


-206 


-311 


B2 


-192 


-275 


CI 


-93 


-131 


Dl 


-162 


-232 


D2 


-108 


-153 



The average half maximum 
beam size in the LR images is 
106 arcsec. 



due to the primary beam. These values for the peak decre- 
ments may be considered limits on the magnitude of any 
decrements in the fields for comparisons with models pre- 
dicting S-Z source counts in blank sky regions. 

It may be noted here that the peak decrements listed in 
Table 5 for the fields A2, Dl and, in particulax, B2 appear 
similar to the decrement listed for field Bl. However, the 
peak decrements seen in the A2, Dl and B2 fields are all 
located almost 3 arcmin from the field centres and appear 
(in images prior to corrections for the attenuation due to the 
primary beam) at levels less than 4 times the thermal rms 
noise. On the other hand, the decrement in the Bl field is 
seen at a level 7.5 times the rms thermal noise and is located 
0.5 arcmin from the field centre. 

The sky area searched for decrements, in Table 5, is six 
fields of size 3 arcmin in radius, corresponding to a total of 
0.05 degrees'^. The limit to the magnitude of decrements in 
this sky area is about 0.2 mjy at 8.7 GHz. 

The Sunyaev-Zeldovich decrement in the CMB bright- 
ness temperature has been detected along the lines of sight 
to several clusters of galaxies (see Birkinshaw 1998 for a 
recent review). For example, an intorferomctric imago of 
the higfi-redshift cluster CL 0016+16 (z = 0.541) with a 
beam of arcmin size showed a decrement corresponding to 
—426 /uK, at lower redshifts, Abell 773 has been seen to show 
a -312 ^iK dip (Carlstrom, Joy & Grego 1996). The S-Z 
decrements detected towards clusters of galaxies in target- 
ted observations have yielded central decrements of several 
hundred /uK. The decrement observed along lines of sight 
through the centre of a cluster depends on the integral of 
the thermal gas pressure along the path and is independent 
of the distance to the cluster; therefore, it is unlikely that 
any rich cluster of galaxies lies in the surveyed region. 

Several authors have estimated the CMB anisotropics 
caused by cosmological populations of clusters of galaxies 
(Markevitch et al. 1994; Bartlett & Silk 1994; Barbosa et al. 
1996). Bartlett et al. (1998) estimate the S-Z source counts 
at 0.75 mm wavelength adopting the Press-Schechter for- 
malism for the mass function and normalizing the matter 
density power spectrum to the present day cluster abun- 
dances as measured by the X-ray temperature distribution 
function. The limit from the ATCA observations translates 
to a limit of 30 mJy flux density at 0.75 mm wavelength. In 
an fio = 1 universe, we may expect about 0.1 sources in the 
ATCA survey region; in an Slo = 0.2 universe the expecta- 
tion is about 1 source. The ATCA survey places constraints 
on models with Oo < 0.2. 



Two separate interferometric images made of sky re- 
gions devoid of obvious clusters appear to have detected 
decrements. Jones et al. (1997) report a —410 jxiy integrated 
decrement in a Ryle-telescope (RT) imago, containing a pair 
of quasars, made with a 110" x 175" beam at a frequency of 
15 GHz. A second field showed no decrement. These corre- 
spond to one ~ 20 mJy S-Z source at 0.75 mm wavelength 
in 0.034 degree^ sky area. Richards et al. (1997) report a 
—27 ^Jy integrated flux density feature in their 30" resolu- 
tion Very-large-array (VLA) images at 8.4 GHz. A separate 
VLA field observed previously showed no such decrement. 
The VLA results correspond to one 4 mJy source at 0.75 mm 
in a sky area of 0.018 degree^. These limits, as well as re- 
sults from the SuZIE (Church et al. 1997) and and OVRO 
RING experiments (Myers et al. 1993) are compared with 
the expected S-Z source counts by Bartlett et al. (1998). 

The constraints imposed on the source counts by the 
ATCA survey are the same as those inferred from the other 
experiments. It may be noted, however, that (i) the ATCA 
survey region considerably exceeds the sky areas covered 
in the other interferometric imaging surveys, and (ii) the 
ATCA survey is truly a 'blank field' search as against the 
RT images which were made towards fields known to have 
quasars. We would also like to point out that the residual 
variance in the fields are dominated by, and consistent with, 
the expectations from residual unsubtracted confusion: con- 
straining the S-Z source counts to be limited to the allowed 
residual variance, after accounting for the variance from 
weak unsubtracted confusion, will be expected to provide 
better limits on S-Z sources and hence on allowed cosmolog- 
ical models. 

6 GAUSSIAN-FORM CMB ANISOTROPIES 

Some CMB anisotropics may be modelled as being Gaus- 
sian random in their temperature distribution. As discussed 
earlier, these anisotropics are usually (and completely) de- 
scribed by the power spectrum which is the distribution of 
l{l + l)Ci in i-spacc. The complex spectral amplitudes at 
different multipole orders are assumed to be random phase 
and the fluctuations on different scales are assumed to be 
uncorrelated. 

Assuming further that the anisotropy spectrum is flat 

implies that l{l + l)Ci is a constant in Z-space. The spectrum 
is then completely specified by its normalization which may 
be parametrized as 

l{l + m = (^){Qnatf. (4) 

The flat-spectrum CMB anisotropy results in a sky variance 

given by 

The normalization Qfiat is referred to as the quadrupole 
normalization because (AT)^ = {Qfiatf for / = 2 implying 
that Qfiat is the rms of the quadrupole anisotropy. 

6.1 ATCA telescope filter function 

In any CMB anisotropy experiment, the instrument 
characteristics, observing strategy and data analysis method 



ATCA survey for CMB anisotropies 13 



fen 




2000 



4000 



6000 



8000 



Multipole order I 



Figure 6. ATCA telescope filter function corresponding to tfie 
sky variance estimate tliat was derived as a weighted mean of the 
variances in the six LR residual images. 



together determine a 'telescope filter function' (TFF) which 
is an /-space filter Fi that defines the sensitivity of the ob- 
servation to the CMB anisotropy power in different multi- 
pole orders. The ATCA LR images at 8.7 GHz have a beam 
of about 2 arcmin and this scale corresponds to multipole 
orders of several thousand. In the large-Z limit, the LR im- 
age rms AS expected from a flat-spectrum CMB anisotropy 
spectrum may be approximately written in the form 



(6) 



where Fi has units (Jy beam~^ K~^)'^. As discussed in Ap- 
pendix A of Subrahmanyan et al. (1998), the TFFs vary 
across the synthesized image and may be computed at any 
image pixel position {6o, <j)o) from the 'product beam' b{9, (f>), 
which is the product of the primary beam centred at the an- 
tenna pointing centre and the synthesized beam centred at 
{do,<i)o). The TFF is given by 



(7) 



where 



{Olcos£,-\-<l>lsin£,) 



(8) 



The sky variance was estimated as a weighted average of 
the variances in the six LR residual images; the weights used 
wore inversely proportional to the thermal-noise variances. 
Therefore, the synthesized beam used in the computation of 
the TFF was obtained by concatenating all the LR data and 
constructing a synthesized beam pattern by weighting the 
data based on the system temperatures recorded along with 
the individual visibility measurements. Because the rms in 
the individual images were computed by weighting the pixels 
using the primary-beam pattern, the TFFs computed at dif- 
ferent image pixel positions were averaged using the squared 
primary-beam pattern as weights. The resulting ATCA TFF 
is shown in Fig. 6. The filter function Fi peaks at I = 4900 



and has half-maximum values at I = 3520 and 6200. Sum- 
ming over the filter values yields 

, _i \ 2 

/iJy beam 



I 



(9) 



A flat-spectrum CMB anisotropy spectrum with normaliza- 
tion Qfiat would, consequently, be expected to result in an 

image rms 



AS* = 0.844 Qfi. 



(10) 



where Qfiat has units of yuK and AS is in /jjy beam"'^. 

It may be noted here that the aperture illumination 
of the ATCA antenna — apart from the visibility-domain 
(u, ti)-track of the 30.6 m baseline — determines the shape 
of the TFF profile. We have modelled the primary-beam 
pattern in the aperture domain taking into account the cen- 
tral blockage of the Cassegrain optics and a radial aperture 
illumination profile. The parameters were determined from 
measurements of the beam pattern. 

6.2 Polarized CMB anisotropies 

The visibility data on the six fields have been used to 
construct LR images in Stokes Q and U in a search for po- 
larized anisotropies on arcmin scales. The image rms were 
computed separately for the six fields after weighting the 
image intensities by the antenna primary-beam pattern and 
these values are listed in Table 6. As discussed earlier, the 
expected image rms due to thermal noise may be derived 
from the corresponding Stokes V images and these values 
are also listed in the table for comparison. 

We have used the thermal-noise variances (variance in 
the corresponding Stokes V images) as weighting factors and 
computed the mean variance in the Stokes Q images to be 
(22.2 /ijy beam~^)^ and that in the Stokes U images to be 
(21.8 ^iJy beam"^)^ The weielited-mcan thermal noise is 
24.1 /LtJy beam^^. There does not appear to be any 'excess' 
variance in the polarized intensity images. 

Examining the images without any correction for the 
primary-beam attenuation, the image intensities within 
the main lobes of the primary beams (within 7 arcmin 
of the image centres) are seen to be in the range -|-86- 
—81 /iJy beam~^. Simulations of LR images that contain 
exclusively a thermal-noise variance of 24.1 /ijy beam~^ re- 
veals that the image intensities may be expected to span 
the range ±82 /xJy beam~^ in any random sample contain- 
ing 12 images. We conclude that the Stokes Q and U images 
of the six fields do not contain any features apart from those 
expected due to thermal noise. 

Following the procedure described in Subrahmanyan et 
al. (1998), we have simulated sky images assuming that the 
sky is devoid of sources — this is probably a reasonable as- 
sumption for the Stokes Q and U images — and contains 
only fiat-spectrum CMB variance. The polarized intensity 
LR images may be assumed to have contributions from ther- 
mal noise and any polarized CMB fluctuations: because the 
fiuctuations owing to these two contributors would be un- 
correlated, we have separately simulated images containing 
CMB fiuctuations and thermal noise. Variance estimates are 
derived from the simulations — separately for the thermal 
noise and CMB contributions — in the same fashion as for 
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Figure 7. Likelihood functions — derived from simulations — that refer to the case of the Stokes Q and U images where there is no 
significant discrete source confusion. The dotted line represents the probability distribution of the thermal-noise contribution to the 
observed image variance. The dashed line is the distribution expected for the CMB variance contribution assuming Q fi^f = 16 fXK. The 
continuous line is the convolution of these two curves and represents the likelihood of obtaining any observed variance conditional on 
Qflat = 16 mK. 



Table 6. R,ms pixel intensities in Stokes Q, 
U and V images. 



Field 


Stokes Q 


Stokes U 


Stokes V 


A2 


25.2 


28.6 


23.9 


Bl 


24.5 


20.5 


26.9 


B2 


25.4 


20.6 


23.0 


CI 


24.0 


24.0 


25.5 


Dl 


28.5 


15.7 


23.0 


D2 


18.4 


19.8 


22.8 



All rms noise values are in jxjy beam~^ . The 
rms in these Stokes Q and U images were 
computed from pixel intensities weighted by 
the primary beam. 

the observations: primary-beam weighted rms is derived for 
each simulated image and the average variance is computed 
for sets of six images. In Fig. 7 we show the probability 
distributions for the variance contributions as inferred from 
the simulations; the CMB anisotropy has been assumed to 
have Qflat = 16 jjK and the thermal noise in the LR images 
is assumed to have an rms of 24.1 ^Jy beam~^. It may be 
noted here that these distribution functions do not follow 
a simple chi-square functional form because the image vari- 
ances have boon estimated as a weighted sum of the squares 
of the pixel intensities. The mean variance in the sinmlatcd 
CMB-anisotropy images is (13.6 /xJy beam^^)'^ and this is 
consistent with that expected from the derived ATCA TFF 
(equation 10). The convolution of the two distribution func- 
tions one for exclusively CMB anisotropy and the other 
corresponding to purely thermal noise — gives the distri- 
bution expected for the observed images and this is also 



shown in Fig. 7. The result of the convolution represents the 
likelihood of observing any variance value assuming CMB 
anisotropy with Qfiat ~ 16 /iK. 

We use the likelihood-ratio test to derive upper limits on 
any polarized anisotropy in the CMB on arcmin scales. The 
null hypothesis Hi is that Q fiat is zero in the /-space window 
covered by the ATCA TFF and that the observed LR images 
in polarized intensity have only thermal-noise variance. The 
likelihood function 'iP{a^^^ | 0) is the probability that the 
observed variance is ali,^ conditional on Qfiat = and is 
given by the dotted line in Fig. 7. If we make the hypothesis 
IH2 that the CMB anisotropy has a flat-band power Qfiat 
in addition to the expected thermal noise, the likelihood 
function JP{a^i,g \ Qfiat) gives the probability distribution 
of the observed variance. For the particular case of Q fiat = 
16 fjK, this likelihood function is given by the continuous 
curve in Fig. 7. 

The decision variance is the observed variance cr^bs- K 
we make the hypothesis IH2 that the CMB anisotropy has a 
flat-band power exceeding Qfiat, the probability of rejecting 
this hypothesis when it is true is given by the integral of 
the likelihood function JP {agi,^ | Qfiat) up to the decision 
variance CTotg. This integral is referred to as the size a of the 
test of hypothesis 1H2; 1 — a is the confidence with which the 
observations with variance (j^jg reject the hypothesis. The 
probability of accepting the hypothesis IH2 when the null 
hypothesis Hi is true defines the quantity 1 — /3, where /3 is 
referred to as the power of the test. 

The observed variance in Stokes Q and U arc less than 
the expected variance from thermal noise. If we adopt the 
values of the observed variances in Stokes Q and U as the 
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Figure 8. Likelihood functions — derived from simulations which include discrete-source confusion as well as thermal noise — that 
refer to the case of the Stokes I images. The dotted line represents the probability distribution of the combined contribution to the 
image variance from thermal noise + residual confusion. The dashed line is the distribution expected for the CMB variance contribution 
assuming Qfiat = 28 /iK. The continuous line is the convolution of these two curves and represents the likelihood of obtaining any 
observed variance conditional on Q fi^i = 28 //K. 



decision variance to reject hypotheses, the tests will have 
powers 13 < 0.5. If choose to adopt the thermal-noise vari- 
ance of (24.1 fijy beam~^)^ as the decision variance; this 
will yield tests with power f3 = 0.5. With this choice, we 
reject the hypothesis that Q fiat exceeds 16 /nK with 95 per 
cent confidence. Adopting, instead, the observed Stokes Q 
variance of (22.2 /xJy beam~^)^ as the decision variance, we 
conclude that Qfiat < 11 fjK with 95 per cent confidence; 
this test has a power /3 = 0.23. The corresponding limit 
based on the observed variance of (21.8 /iJy beam-i)2 in 
Stokes U is Q flat < 10 ijK with 95 per cent confidence and 
power P = 0.17. 

6.3 Limits on flat-band total-intensity anisotropy 

We derive limits on flat-band CMB anisotropies in total in- 
tensity from the LR images of the six surveyed fields. The 
residual images, shown in Fig. 3, are expected to have contri- 
butions from (i) the telescope thermal noise, (ii) fluctuations 
due to residual unsubtracted weak sources in the fields and 
(iii) any CMB anisotropy. The weighted mean of the vari- 
ances in the six residual LR Stokes I images was derived 
to be (52.0 /uJy beam"'^)^. As discussed earlier, the tele- 
scope thermal noise variance contribution is estimated to 
be (24.1 /xJy beam^^)^ from the LR Stokes V images and, 
therefore, there is a significant variance in the observed im- 
ages apart from the thermal noise. 

The simulations of discrete-source confusion, discussed 
in section 3.2, indicated that if we adopt the nominal source 
counts given in equation 2, the residual images may be ex- 
pected to have a variance of (59 ± 8 /uJy beam~^)^ owing 



to thermal noise and the unsubtracted confusion. This is 
consistent with the observed variance in the residual images 
and we conclude that there is no evidence in the LR images 
for any CMB anisotropy contribution. 

It may also be noted that the HR residual images 
also have a significant variance in excess of the expected 
thermal noise — (25.0 jiJy beam~^)^ as compared to the 
(19.4 fjbJy beam~^)^ thermal noise contribution — and this 
is also consistent with the simulations which predict a vari- 
ance (24.4±1.2 /ijy beam"^)'^ for the combined contribution 
from thermal noise and residual confusion. 

We assume that the sky fluctuations due to any CMB 
anisotropy is uncorrelated with the fluctuations due to ther- 
mal noise and residual confusion. The simulations described 
in section 3.2, which provided random samples of images 
containing the expected contributions from thermal noise 
and residual confusion, were used to derive the expected 
distribution of image variance conditional on Qfiat = 0. 
This distribution, ]P{a'^ijg \ 0), is shown as the dotted line 
in Fig. 8 and represents the likelihood function for the null 
hypothesis Hi . The distribution function for the image vari- 
ance contribution from a flat-band CMB anisotropy with 
Qfiat = 28 /jK is shown in Fig. 8 as a dashed line. The 
convolution of the two yields the continuous line in the fig- 
ure and represents the likelihood function FifT^ts I Qfiat) 
for the particular assumption that Qfiat ~ 28 piK. The ob- 
served variance in the Stokes I images, (52.0 /iJy beam^^)^, 
rejects the hypothesis that Q fiat exceeds 28 /xK with 95 per 
cent confidence. 

The LR images, before any confusion subtraction, have 
an rms that exceeds the thermal noise rms by a factor almost 
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5. After subtracting the confusion models, the residual LR 
images are still confusion limited with an rms that is a factor 
about 2 higher than the thermal noise. Therefore, the esti- 
mation of any CMB anisotropy contribution in the residual 
images is critically dependent on our knowledge of, and un- 
certainties associated with, the source counts at the 10-/xJy 
level. Overestimations of the source-count normalization will 
result in spuriously low limits on any CMB anisotropy. How- 
ever, we do not believe we have overestimated the expected 
discrete source counts because (i) the variance in the resid- 
ual HR images exceeds the expected variance by 0.5a and 
(ii) the source counts derived by Partridge et al. (1997) in 
independent observations nominally exceeds the counts we 
have adopted here at the 10 — /iJy flux density level. 

We have incorporated the uncertainties in the normal- 
ization and index of the differential source counts by sim- 
ulating fields with Gaussian random choices of parameters 
for the source counts; we adopt 1-standard-deviation spreads 
for the parameters corresponding to the uncertainties quoted 
in equation 2. These simulations show, as may be expected, 
a 12 per cent increase in the standard deviation of image 
variances; however, the mean imago variance is biased to a 
larger value: (60.5 /iJy bcam^^)"^. Using this distribution as 
the likelihood function P(CT^5g | 0) of the null hypothesis, 
and adopting the observed variance of (52.0 /iJy beam-i)2 
as the decision variance, we find that the observations rule 
out flat-band CMB anisotropy with Qfut > 25 /^K with 95 
per cent confidence. 

6.4 Discussion 

The 4-year COBE observations of CMB anisotropy (Ben- 
nett et al. 1996) correspond to a flat-band spectral power 
of Qfiat = 18 /iK in the multipole range 2^1^ 20. Sev- 
eral experiments measuring anisotropies in /-space windows 
with centres located at multipoles up to a few hundred have 
reported detections corresponding to much greater spectral 
powers (see the website 

www. sns .ias.edu/ ~max /cmb / experiments . html 
maintained by Max Tegmark for a compilation of experimen- 
tal results). For example, the recent Python V results (Coble 

et al. 1999) report band powers rising from Qfut = 15 to 
63 /jK over the multipole range / = 50 to 200. 

At multipoles exceeding I ~ 10^, there are no reports to 
date of the detection of Gaussian anisotropies in the CMB; 
observations have, however, detected isolated 'holes' in the 
CMB towards clusters of galaxies. As discussed earlier, the 
tentative detections of 'holes' in 'blank fields' are also pre- 
sumably owing to the S-Z effect in distant clusters. The 
SuZIE experiment (Ganga et al. 1997) limit Q fiat < 27 /iK 
at the 2-(7 level in the multipole window / = 1330-3670. 
Rylo telescope observations of the Lynx 2 field (Jones 1998) 
limit Qfiat < 12 pK, or AT/T < 6.6 x 10"®, in the 
window 4300 < I < 8300 with 67 per cent confidence. 
The OVRO NCP experiment (Readhead et al. 1989) limit 
AT/T < 1.7 X 10"'' at an angular scale of 2 arcmin. VLA 
images at 8.7 GHz limit AT/T < 2.0 x 10"''' at an angular 
scale of 1 arcmin (Partridge ct al. 1997). The VLA images 
also limit the sky anisotropy in linearly polarized intensity 
to AT/T < 1.0 X lO""' on the 1-arcmin scale. 

The ATCA observations reported here cover six in- 
dependent regions with a combined sky area of about 



0.05 degrees^. The images limit the CMB anisotropy in the 
multipole range / = 3520-6200 to Qfiat < 25 /iK in total in- 
tensity and 10 /uK in linearly polarized fluctuations. These 
correspond to upper limits of {AT/T) < 6 x 10~® in linearly 
polarized intensity and < 1.4 x 10~^ in total intensity for 
fluctuations with an angular scale of 2 arcmin. The limits 
on polarized intensity fluctuations are the best to date; the 
limits on total intensity fluctuations are comparable to those 
set by other experiments on arcmin scale. These results are 
an improvement over the results published in Subrahmanyan 
ct al. (1998) in that the inevitable confusion variance contri- 
bution to arcmin-scale sky fluctuations at 8.7 GHz has now 
been accounted for. The small-angle anisotropy experiments 
have consistently set limits that clearly establish the down- 
turn in the CMB anisotropy spectrum between I ~ 10^ and 
10^ 



7 CONCLUSIONS 

We have made sensitive observations of six 'blank' fields 
with the aim of examining the arcmin scale anisotropy in the 
CMB. After subtracting confusion models derived by higher 
resolution data, the residual images are consistent with the 
expected thermal noise and residual confusion from weak 
undetected sources. We do not believe we have detected any 
S-Z decrements in the blank fields. This database is useful for 
deriving constraints (i) on the gas in cosmological distribu- 
tions of clusters of galaxies, and (ii) on the CMB anisotropy 
power spectrum in total intensity and polarization at arcmin 
scales. 

The variance in the images at 8.7 GHz is consistent with 
the expectations from thermal noise and the source counts 
at /ttJy flux density levels derived by Windhorst et al. (1993). 

Assuming Gaussian- form CMB anisotropy with a 'flat' 
spectrum, we derive 95 per cent confidence upper limits of 
Qfiat < 10 11 /iK in polarized intensity and Qfiat < 25 /iK 
in total intensity. These correspond to upper limits on 
(AT/Tcmb) of 5.7 6.2 x 10"® in polarized intensity and 
1.4 X 10~® in total intensity, where Tcmb = 2.73 K and 
AT = [/(/ + l)Ci/27r]5. The ATCA filter function peaks at 
/ = 4700 and has half maximum values at I = 3350 and 
6050. 

The ATCA observations are confusion limited and im- 
provements in detection capability require observations at 
higher frequencies where the confusion is less. The ATCA 
is currently being upgraded to operate in the 16-25 GHz 
band and a north spur is being added to the railtrack to 
allow 2-dimensional close-packed array configurations. As 
discussed in section 2, we expect the confusion noise to re- 
duce as AS ~ i^"^ '^, and observations at 18 GHz would 
have a confusion rms noise reduced by a factor 5 as com- 
pared to the current observations at 8.7 GHz. The angu- 
lar scale of the imaging would be reduced by a factor 2 to 
about 1 arcmin. High-brightness-sensitive observations with 
an ultra-compact ATCA at 18 GHz may improve limits on 
CMB anisotropy on arcmin scales. In particular, the upgrade 
would be useful in resolving the sub-structure in the X-ray 
gas of higli-redshift clusters — detected in future blank-field 
S-Z surveys — by imaging their S-Z decrement. 

The ATCA CMB anisotropy research website is located 
at www.atnf.csiro.au/research/cmbr/cmbr_atca.html. 
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